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PREFACE 


Pacific-Sierra  is  performing  a  study  of  possible  military  applica¬ 
tions  for  aii  electromagnetic  waveguide  in  the  lithosphere.  This  report 
presents  the  results  of  the  first  phase  of  that  study;  viz.,  a  determina¬ 
tion  of  the  data  rates  and  transmission  ranges  that,  based  on  currently 
used  models  of  the  lithosphere,  might  be  achievable  with  reasonable  power 
expenditure  and  antenna  burial  depths.  These  results  form  the  scientific 
basis  for  the  second  phase  of  the  study,  which  is  to  determine  those 
missions  for  which  use  of  a  lithospheric  communication  link  could  be 
feasible . 

The  motivation  for  this  work  stems  from  the  integrated  studies  of 
crustal  properties  that  have  been  conducted  for  a  number  of  years  by 
the  Earth  Physics  Program  of  the  Office  of  Naval  Research.  This  effort 
was  brought  into  focus  by  n  Symposium,  jointly  sponsored  by  ONR  and  CIRES, 
held  in  July  1970.  The  proceedings  of  this  symposium  were  published  in 
AGU  Monograph  //14,  The  structure  and  Fhysioal  Properties  of  the  Earth’s 
Crust  (2971).  Subsequently,  two  workshops,  funded  by  ARPA  and  organized 
by  ONR,  were  held  at  the  Colorado  School  of  Mines,  in  February  1972,  to 
define  an  experimental  program  to  ascertain  the  existence  tf  the  litho¬ 
spheric  waveguide. 
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SUMMARY  AND  CONCLUSIONS 


This  report  presents  results  of  an  analysis  of  the  capabilities 
that  might  he  achieved  with  a  lithospheric  communication  system.  The 
analysis  uses  several  available  models  of  continental  lithospheric 
conductivity  as  inputs  to  calculations  of  transmission  characteristics 
and  possible  communication  system  parameters.  Because  of  prevailing 
uncertainties,  significant  differences  exist  among  the  various  models; 
and  a  wide  range  of  future  system  capabilities  is  found  to  lie  within 
the  realm  of  possibility.  Full-wave  modal  solutions  that  take  full 
account  of  vertical  conductivity  gradients  are  used  to  calculate  attenua¬ 
tion  rates,  field-strength  depth-profiles,  and  transfer  functions 
relating  the  noise  fields  at  depth  to  those  measured  at  the  surface. 

Also,  estimated  values  of  required  system  power  are  given  as  a  function 
of  data  rate,  transmission  range,  and  receiver  burial  depth.  Carrier 
frequencies  from  100  Hz  to  100  kHz  are  considered.  Ic  could  well  prove 
worthwhile,  however,  to  extend  the  analysis  to  frequencies  outside  of 
the  ELF/VLF/LF  bands  treated  herein. 

For  all  models  used,  the  calculated  attenuation  rates  are  con¬ 
siderably  larger  than  for  above-ground  transmission,  which  depends  on 
propagation  in  the  earth- ionosphere  waveguide.  In  spite  of  this  fact, 
reasonably  large  transmission  ranges  might  be  possible  in  the  lithosphere 
because  the  atmospheric  noise  fields  at  depth  are  very  small,  having 
suffered  heavy  attenuation  in  propagating  downward  from  the  earth's 
surface . 

The  major  factor  affecting  system  feasibility  is  whether  water- 
saturated  microcracks  exist  in  all  major  rock  types  at  depths  greater 
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than,  say,  5-to-8  km.  If  such  microcracks  exist  (as  assumed  in  one 

-5  -6  , 

model),  and  minimum  conductivities  in  the  crust  exceed  10  -  10  mhos/m, 

then  the  calculations  show  clearly  that  a  practical  communication  system 
is  not  feasible.  For  this  case,  the  attenuation  rates  are  so  high  that 
only  very  short  transmission  ranges  could  be  achieved,  even  if  large 
amounts  of  power  were  expended.  Field  experiments  will  be  needed  to 
determine  with  certainty  whether  saturated  microcracks  are,  in  fact, 
present  at  the  depths  of  Interest.  However,  valid  reasons  exist  for 
optimism  that  fluids  (or  fluid  films)  should  not  be  an  Important  factor 
at  depths  greater  than  several  kilometers,  and  that  conductivity  profiles 
(at  depth)  based  on  laboratory  data  for  dry  rocks  more  closely  approxi¬ 
mate  actual  conditions.  Further,  recent  laboratory  data  indicate  that 
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the  conductivity  of  dry  rock  in  the  crust  could  be  much  lower  (10  to 
10  ^  mhos/m)  than  previously  believed  (10  mhos/m).  It  is,  of  course, 
difficult  to  relate  data  obtained  from  laboratory  samples  to  conductivities 
of  rock  in  situ.  Nonetheless,  calculations  using  conceptual  profiles 
based  on  the  assumption  of  dry  basement  rock  yield  results  that  are 
quite  encouraging.  For  example,  for  data  rates  of  a  few  tens  of  bits- 
per-second  and  power  expenditures  of  1-to-lO  megawatts,  transmission 
ranges  of  from  about  1000  to  many  thousands  of  kilometers  (depending 
on  the  profile  used)  appear  to  be  a  possibility. 

Tlie  calculations  show  clearly  that  much  larger  transmission  ranges 
can  be  achieved  in  the  LF  band  (30-100  kHz)  than  at  the  lower  fre¬ 
quencies  considered.  This  effect  is  due  mainly  to  the  heavy  attenua¬ 
tion  suffered  by  downward  propagating  atmospheric  noise  at  the  higlier 
The  calculations  also  show  that  the  maximum  transmission 


frequencies . 
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ranges  are  somewhat  sensitive  to  the  details  of  the  model  profiles. 

For  example,  for  two  profiles  exhibiting  the  same  minimum  conductivity, 
tlie  transmission  range  (for  a  specified  power,  etc.)  can  vary  by  a 
factor  of  two  depending  on  the  conductivity  deptli-gradients ,  and  tlie 
precise  jocation  of  any  sharp  transition  layers  that  miglit  exist. 

The  depth  to  which  a  receiver  must  be  buried  to  achieve  satisfac¬ 
tory  performance  is  a  particularly  important  parameter,  since  borehole 
d’'illing  costs  could  be  a  large  fraction  of  tlie  total  system  cost. 

The  transmission  ranges  given  above  apply  when  the  receiving  antenna 
is  'ccated  near  the  "center"  of  the  waveguide — assumed  to  be  about 
ten  kilometers  deep  for  the  conductivity  profiles  used  in  this 
report.  For  profiles  representative  or  regions  stripped  of  highly 
conduct! /e  sedimentary  layers,  the  calculations  show  that  receiver 
depths  as  shallow  as  3  or  4  km  could  be  used  at  only  a  15-co-20  percent 
penalty  in  transmission  range.  One  implication  of  this  fact  is  that — 
previous  negative  results  notwithstanding — a  meaningful  propagation 
experiment  could  be  carried  out  with  relatively  shallow  boreholes, 
provided  that  the  region  was  carefully  selected.  Frequencies  of 
several  tens  of  kilohertz  would  be  preferable,  and  care  sliould  be 
taken  to  insure  a  reasonable  transmitter  efficiency.  Of  course, 
deeper  borehcles  would  be  needed  in  regions  having  thick,  highly 
conducting  overburdens. 

Another  type  of  propagation  experiment,  not  involving  a  trans¬ 
mitter,  has  appeal  for  ascertaining  the  existence,  or  non-existence, 
of  a  lithospheric  waveguide.  There  are  regions  of  the  sea  floor  where 
little  or  no  sedimentation  exists,  and  where  energy  propagating  in  a 
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lithospheric  waveguide  could  conceivably  approach  the  ocean's  bottom 
from  below.  In  deep  water,  the  preferred  propagation  path  (to  the 
bottom)  for  atmospheric  noise  would  be  via  a  lithospheric  duct,  rather 
than  downward  from  the  ocean’s  surface,  due  to  the  opacity  of  sea  water 
to  all  frequencies  higher  than.  say.  a  few  Hertz.  Consider  an  experi¬ 
ment  where  atmospheric  noise  was  measured  as  a  function  of  the  depth 
(in  the  water)  of  the  receiving  antenna.  For  shallow  receiver  depths, 
this  noise  would  decrease  as  the  receiver  depth  was  Increased.  If. 
however,  the  noise  (e.g..  Schumann  resonances)  was  then  found  to 
increase  as  the  receiver  apMoached  or  contacted  the  bottom,  consider¬ 
able  credence  would  be  given  to  the  presumed  existence  of  a  crustal 

waveguide . 
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REFERENCES 


I.  INTRODUCTION 


More  than  a  decade  has  passed  since  the  existence  of  an  electro¬ 
magnetic  waveguide  in  the  earth's  crust  was  first  postulated  (e>g>i 
Wait,  19[>4;  IWheeler,  !961).  The  assumption  was  that  conductive,  wet 
surface  layers  and  conductive,  hot  mantle  layers  enclose  a  resistive 
zone  of  dry  basement  rocks.  As  originally  conceived,  the  resistive 
zone  occupied  depths  between  roughly  10  and  30  km.  Interest  in  ascer¬ 
taining  the  existence,  or  non-existence,  of  this  lithospheric  waveguide 
stemmed  from  its  potential  use  as  a  communication  channel.  In  the 
intervening  years  enthusiasm  about  the  prospect  of  a  crustal  waveguide 
waned,  partially  because  of  the  apparently  negative  outcome  of  propaga¬ 
tion  experiments  carried  out  between  1963  and  1966  (  ■•  '  .'  ) . 

In  retrospect,  it  appears  that,  due  to  experimental  limitations,  the 
data  thus  obtained  could  not  have  either  confirmed  or  denied  the  exis¬ 
tence  of  a  waveguide  in  the  lithosphere. 

Little  is  known  in  detail  about  the  precise  nature  of  the  materials 

and  physical  properties  of  the  earth's  crust  at  depths  where  the  resistive 

basement  rocks  are  presumed  to  exist.  The  available  data  base  is  by 

no  means  complete,  and  additional  field  measurements  are  needed  if  firm 

values  are  to  be  assigned  to  key  parameters.  However,  reinterpretation 

of  available  field  data  has  indicated  that  conductivities  of  10  mhos/ra 

or  lower  could  exist.  Moreover,  recent  laboratory  experiments  on  the 

conductivity  of  dry  rock,  and  recent  evidence  based  on  minerological 

considerations,  suggest  cause  for  optimism  that  conductivities  lower 
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than  previously  believed  might  exist  in  the  crust  (e.g.,  '  10  mhos/m). 
Theoretical  analyses  show  that  the  corresponding  attenuation  rates  for 


electromagnetic  waves  could  be  low  enough  to  permit  communication 
over  useful  distances.  For  these  reasons,  there  has  been  renewed 
interest  in  determining  the  electrical  properties  of  the  lithosphere 
and  in  assessing  its  potential  utility  as  a  communication  channel. 
Workshops  have  been  held,  and  recommendations  made  that  interdisci¬ 
plinary  field-measurement  programs  be  undertaken  {Keller,  1971;  Haler,, 
'J72) . 

This  report  presents  results  of  an  analysis  of  capabilities  that 
might  be  achieved  with  a  lithospheric  communication  system.  As  might 
be  expected  in  view  of  prevailing  uncertainties,  a  wide  range  of  future 
system  capabilities  is  found  to  lie  within  the  realm  of  possibility. 
However,  the  results  give  guidance  as  to  whether  it  is  reasonable  to 
expect,  on  the  basis  of  current  knowledge,  that  a  lithospheric  wave¬ 
guide  would  provide  the  basis  for  useful  communication  links. 

The  lithospheric  conductivity  models  used  in  the  calculations  are 
given  and  discussed  in  Sec.  II,  as  are  models  of  atmospheric  and  therma 
noise.  The  mathematical  methods  used  are  described  in  Sec.  Ill,  and 
numerical  results,  obtained  from  Pacific-Sierra's  long-wave  propagation 
code  are  given  in  Sec.  IV. 


II.  MODELS  OF  ELECTRICAL  CONDUCTIVITY  AND  NOISE  IN  THE  LITHOSPHERE 

The  electromagnetic  transmission  properties  of  the  lithosphere 
must  be  calculated  in  order  to  assess  its  utility  as  a  communication 
channel.  This  calculation  requires  a  knowledge  of  the  depth-profiles 
of  the  electrical  properties  (primarily  conductivity)  of  the  earth's 
crust  and  upper  mantle,  as  well  as  a  knowledge  of  the  degree  of  lateral 
homoegneity  of  these  properties.  In  addition,  tl  noise,  environment 
at  the  receiver  must  be  known  to  determine  the  relationships  among 
radiated  power,  transmission  range,  and  data  rate.  Accordingly,  the 
models  of  conductivity  and  noise  spectrum  used  as  inputs  to  the  calcula¬ 
tions  of  Sacs.  Ill  and  IV  are  discussed  below. 

CONDUCTIVITY  DEPTH-PROFILES 

Considerable  uncertainty  exists  as  to  the  conductivity  of  the 
earth  at  depths  greater  than  several  kilometers.  This  uncertainty  is 
due  to  practical  difficulties  encountered  in  using  surface-based  sound¬ 
ing  techniques  to  detect  low  conductivities  at  depth,  and  an  absence 
of  sufficiently  deep  boreholes  for  direct  measurements.  However, 
nominal  models  for  conductivity  depth-profiles  in  the  litho-phere 
have  been  suggested  by  various  geophysicists.  As  might  be  expected 
in  view  of  the  prevailing  uncertainties,  significant  differences  exist 
among  currently  available  model  profiles.  This  state  of  affairs  is 
unfortunate,  since  the  results  given  in  Sec.  IV  show  that  the  trans¬ 
mission  properties  of  the  lithosphere  depend  strongly  on  the  minimum 
value  of  conductivity  in  the  crust,  as  well  as  on  the  depth-gradients 


of  the  conductivity. 
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It  is  necessarily  beyond  the  scope  of  this  study  to  determine  the 
correct  depth-profile  of  conductivity  In  the  lithosphere.  Indeed,  the 
main  goal  of  the  Crustal  Studies  norkshop  (Hales,  et  al.,  1972)  was  to 
consider  and  suggest  experimental  procedures  by  which  profiles  could  be 
determined.  The  approach  taken  In  this  report  Is  to  use  currently 
available  conductivity  oflles,  which  In  our  opinion  are  reasonable 
In  the  context  of  present  knowledge,  as  Inputs  to  our  calculations. 

This  approach  provides  the  most  reliable  estimates  that  can  now  be  made 
of  the  expected  performance  of  a  lithospheric  communication  system,  and 
also  indicates  the  sensitivity  of  lithospheric  transmission  properties 
to  variations  in  the  profiles.  Guidance  Is  thereby  obtained  as  to 
whether  more  extensive  experimental  determinations  of  lithospheric 
conductivities  are  likely  to  reveal  the  existence  of  a  useful  crustal 
waveguide . 

For  a  discussion  of  evidence  for  (and  against)  the  existence  of 
an  electromagnetic  waveguide  In  the  lithosphere,  the  reader  Is  referred 
to  AGV  Monograph  14  (Heaaock,  ed. ,  1971)  and  Crustal  Studies  Workshop 
Report  (Hales,  1972).  Briefly,  the  conductivity  of  rocks  near  the 
earth's  surface  Is  relatively  high  due  to  the  presence  of  water  In 
pores  and  cracks.  As  the  depth  Is  Increased,  the  contribution  of  pores 
and  cracks  to  the  conductivity  Is  diminished  due  to  Increasing  pressure 
the  conductivity  will  therefore  decrease  with  Increasing  depth.  There 
Is,  however,  a  competing  effect;  viz.,  the  temperature  Increases  by 
some  10'’-to-30‘’K  for  each  additional  kilometer  of  depth.  Ultimately, 
conductivity  due  to  thermal-activated  carriers  In  hot  rock  should  be¬ 
come  dominant.  The  conductivity  should  then  Increase  with  Increasing 
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depth.  There  is  thus  reason  to  expect  a  layer  of  relatively  low 
conductivity  to  exist  between  a  surface  layer  (overburden)  that  is 
highly  conductive  because  it  contains  fluid  and  a  sub-basement  that 
is  highly  conductive  because  it  is  hot.  This  poorly  conducting  layer, 
if  it  exists,  would  comprise  the  lithospheric  waveguide  that  is  the 
subject  of  this  study.  Geophysicists  have  given  much  attention  to 
ascertaining  the  minimum  value  of  conductivity  in  the  waveguide,  since 
this  value  strongly  Influences  the  transmission  properties.  Presently, 
the  value  of  minimum  conductivity  is  uncertain  by  several  orders  of 
magnitude . 

For  the  purpose  of  making  calculations,  the  so-callcd  step- 
function"  representation  of  the  conductivity  profile  is  the  simplest. 
Here,  the  lithosphere  is  assumed  to  form  a  uniform  parallel-plate 
waveguide,  bounded  on  the  top  by  a  highly  conducting  overburden  and  on 
the  bottom  by  a  highly  conducting  sub-basement.  The  region  between 
these  two  sharp  boundaries  contains  rock  of  low  conductivity.  Cpiee 
and  Wait  (1971)  performed  an  extensive  and  useful  parametric  study  of 
the  transmission  properties  of  a  parallel-plate  model  of  the  lithosphere. 
Specifically,  they  calculated  attenuation  rates  and  excitation  functions 
as  functions  of  wave  frequency,  waveguide  width,  waveguide  conductivity, 
and  boundary  conductivity.  Although  the  step-function  representation 
of  the  conductivity  profile  provides  insight  into  tha  dependence  of  the 
transmission  on  several  key  parameters,  we  have  chosen  to  use  more 
detailed  profiles.  This  choice  was  made  because  1)  the  signal  attenua¬ 
tion  depends  on  the  depth-gradients  of  the  conductivity,  which  cannot 
be  included  in  a  step-function  profile;  and  2)  the  depth  to  which 
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antennas  must  be  burled  to  provide  adequate  reception  (or  radiation) 
cannot  be  calculated  satisfactorily  from  a  step-iunction  model.  To 
ascertain  the  depth  dependences  of  the  signal  and  noise  fields,  finite 
conductivity  gradients  must  be  included  in  the  model. 

A  major  source  of  uncertainty  is  the  water  content  of  rocks  at 
depth.s  greater  than  a  few  kilometers.  One  view  holds  that  microfractures 
are  not  present  in  rocks  at  depths  that  have  never  been  mechanically  un¬ 
loaded  from  lithostatic  compression.  Such  rocks  would  be  dry  and  have 
a  very  low  electrical  conductivity.  Another  view  is  that  water  films 

•k 

in  microcracks  must  be  present  at  all  depths  in  the  crust  {Brace,  1971). 
Hypothetical  conductivity  profiles  that  illustrate  the  implications  of 
the  above  two  points  of  view  have  been  prepared  by  levin  (1971)  and  are 
shown  in  Fig.  1.  Each  profile  exhibits  a  "waveguide"  at  depths  between 
about  10  and  30  km.  The  conductivities  at  shallow  depths  are  based  upon 
available,  shallow  geophysical  measurements;  those  below  about  30  km 
are  estimated  from  assumed  temperature  gradients  in  the  rock.  The  profiies 
shown  in  Fig.  1  differ  only  in  the  "waveguide"  region  between  10  and  30  km. 
The  segments  labeled  "wet"  and  "dry"  are  nominal  representations  based 
upon  early  labora-.ory  measurements  of  "wet"  and  "dry"  rocks,  respectively. 
Under  the  hypothesis  of  "wet"  rock,  the  minimum  conductivity  shown  is 
about  10  ^  mhos/m;  and  for  dry  rock,  about  10  ^  mhos/m.  The  profiles 
in  Fig.  1  are  intended  only  to  illustrate  the  range  of  uncertainty  in 
the  expected  minimum  conductivity,  and  are  nominal  rather  than  detailed 

*Whether  such  saturated  microcracks  exist  in  all  major  rock  types 
in  situ  is  still  undetermined.  However,  for  several  reasons,  there 
appears  to  be  increasing  optimism  that  dry,  low-conductivity  zones  can 
exist  in  the  crust  (see,  e.g..  Hales,  2972;  p.  28). 
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representations .  Further,  as  discussed  below,  recent  laboratory  measure- 
ments--not  available  when  Levin’s  profiles  were  prepared— indicate  that 
the  minimum  concuctivity  could  well  be  considerably  less  than  10  ^  mhos/m. 

Of  the  currently  available  conductivity  profiles,  the  one  based 
on  the  most  extensive  measurements  is  probably  that  shown  in  Fig.  2 
{Keller  as  reported  by  GallccDa  and  Haidle,  1972).  Although  uncertain 
in  several  respects,  this  profile  represents  a  compilation  of  a  good 
deal  of  geophysical  and  laboratory  -data.  The  conductivities  shown  for 
relatively  shallow  depths  are  based  upon  well-log  data,  taken  in 
Appalachia,  that  clearly  indicate  a  decreasing  conductivity  down  to  at 
least  several  kilometers.  This  decreasing  trend  in  conductivity  is 
believed  to  be  controlled  by  the  rate  at  which  fractures  close  due  to 
increasing  overburden  pressure.  Below,  say,  3-to-5  km,  the  decrease  in 
conductivity  was  extrapolated,  in  a  manner  consistent  with  available 
field  and  laboratory  data,  to  a  minimum  value  referred  from  recent  field 
measurements.*  Below  about  lO-to-12  km,  the  conductivity  is  assumed  to 
be  dominated  by  high-temperature  mineral  conduction.  The  values  shown 
for  depths  below  lO-to-12  km  in  Fig.  2  are  based  on  laboratory  measure¬ 
ments  of  high-temperature  rocks.  The  assumed  temperature  scale  is  shown 
on  the  right-hand  axis  of  Fig.  2,  and  corresponds  to  a  nominal  thermal 

gradient  of  Ib^/km. 

In  addition  to  temperature.,  the  conductivity  at  depth  also  depends 
on  mineralogy.  Note  that  in  Fig.  2  Keller  has  assumed  a  rapid  transition 

I 

^  *For  a  discussion  of  the  various  methods  of  probing  the  electrical 

properties  of  the  crust,  and  a  discussion  of  the  difficulty  in  obtaining 
and  interpreting  da* a  on  deep  layers  of  low  conductivity,  see  Keller 
^  (1971). 

\ 
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in  conductivity  at  about  14  km.  The  existence  of  such  a  transition  is 

inferred  from  seismic  data  and  certain  electromagnetic  field  measurements. 

The  lA-km  depth  at  which  the  transition  is  assumed  to  occur  is,  of  course, 

nviminal,  and  represents  an  average  value  for  the  United  States.  In  some 

regions,  assumed  transition-layer  depths  as  large  as  20-to-25  km  are 

consistent  with  the  available  data  (Keller ,  197?).  Another  characteristic 

of  the  Keller  profile  is  that  the  conductivity  at  the  earth's  surface 

-4 

is  taken  to  be  about  10  mhos/m — a  very  low  value.  The  use  of  a  higher 

“2  “3 

surface  conductivity  (say,  10  to  10  mhos/m)  would  not  change  the  con¬ 
clusions  of  this  report  significantly,  since  the  transmission  properties 
of  the  waveguide  depend  mainly  on  the  properties  of  the  crust  at  depths 
below  1  or  2  km. 

Recent  laboratory  measuraments  of  the  conductivity  of  dry  rock 
reported  by  Houslcy  (197?)  suggest  that  the  minimum  conductivity  of  the 
crust  might  be  considerably  less  than  illustrated  in  Figs.  1  and  2. 

He  argues  that  most  previous  laboratory  data  tend  to  give  a  misleading 
picture,  of  lithospheric  conductivity  at  depths  greater  than,  say,  10  km. 
Rocks  in  situ  should  have  much  lower  conductivity  than  rocks  that  have 
been  exposed  to  the  atmosphere  and  that  have  had  microcracks  closed  by 
the  application  of  high  pressures  in  the  laboratory.  In  the  latter 
case,  thin  films  of  interstitial  fluid  remain  and  raise  the  conductivity 
by  a  considerable  amount.  In  addition,  Housley  has  found  that  the 
partial  pressure  of  oxygen  in  the  atmosphere  surrounding  the  sample 
must  be  maintained  at  a  level  corresponding  to  that  in  the  crystal 
at  the  temperature  and  pressure  for  the  depth  being  simulated.  The 


data  thus  obtained  indicate  conductivities  orders  of  magnitude  less 
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than  previously  measured  for  dry  rocks  in  the  laboratory .  If  one  uses 

these  new  conda':tivity  data,  and  follows  the  same  reasoning  (see  above) 

that  led  to  the  profile  shown  in  Fig.  2,  the  resultant  profiles  are  as 

shown  in  Fig.  3.  These  profiles  are  identical  to  Keller's  (Fig.  2)  for 

depths  less  than  about  10  km,  since  the  same  data  and  assumptions  are 

-8 

used.  At  greater  depths,  however,  much  lower  conductivities  (<10  mhos/m) 
are  indicated,  due  to  the  new  laboratory  data  described  above.  The 
three  profiles  shown  correspond  to  geotherms  in  three  different  heat- 
flow  provinces,  as  reported  by  BlockwstZ  (1971),  The  thermal  gradients 
used  were  BR  (-24'’/km),  EUS  (~l4“/km),  and  SN  (~9“/km)  as  compared 
with  the  Ib'/km  used  for  the  Keller  profile  (Fig.  2). 

It  must  be  emphasized  that  the  Fig.  3  profiles  (or  any  others) 
are  by  no  means  definitive,  and  much  experimental  work  remains  (some 
of  which  is  now  underway)  before  a  satisfactory  understanding  of  the 
conductivity  of  in  situ  rocks  can  be  attained.  Great  uncertainty  is 
involved  in  using  laboratory  data  to  infer  the  conductivity  of  the 
crust.  Further,  the  data  reported  by  Housley  (197Z)  were  for  only  a 
single  sample  (olivine);  but  lithologies  are  not  constant  in  nature. 
However,  these  data  do  give  encouragement  that  the  minimum  conductivity 
of  the  crust  might  be  lower  than  previously  believed. 

Clearly,  the  above  discussion  indicates  selecting  a  single  "best" 
profile  for  making  calculations  would  be  impossible.  Indeed,  Fig.  3 
alone  shows  the  large  variations  that  would  occur  among  different 
provinces.  The  Levin  profile  (Fig.  1)  exhibits  rather  abrupt  changes 
in  c inductivity  at  depths  of  about  2  km  and  7  km.  The  Keller  profile 

*BR  (Basin  and  Range),  EUS  (Eastern,  US),  SN  (Sierra  Nevada). 
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(Fig.  2)  exhibits  a  near-discontin'- ity  at  around  14  km;  whereas  the 

A 

Housley  profiles  (Fig.  3)  are  quite  smooth  at  all  depths.  Field 
measurements  are  needed  to  determine  which  of  the  differences  between 
these  (and  other)  profiles  are  artificial,  and  which  correctly  account 
for  real  differences  among  various  regions  of  the  earth's  crust.  It 
would  not  be  practical,  or  very  meaningful,  to  undertake  detailed  trans¬ 
mission  calculations  for  all  available  conceptual  profiles.  In  Sec.  IV, 
rather  complete  computational  results  are  given  for  the  Keller  and 
Housley  BR  profiles;  and  sample  results  for  other  profiles,  to  indicate 
important  trends.  When  more  accurate  conductivity  measurements  become 
available,  the  methodology  developed  here  can  be  applied  tn  obtain 
more  reliable  estimates  of  the  transmission  properties. 

ATMOSPHERIC  AND  THERMAL  NOISE 

The  noise  environment  at  a  buried  receiver  consists  of  thermal 
noise,  which  depends  upon  the  temperature  distribution  in  the  crust,  and 
atmospheric  noise.-  which  propagates  downward  from  the  earth's  surface. 

We  use  the  temperature  profile  given  in  Fig.  2  to  compute  thermal  noise. 
As  discussed  in  Sec.  Ill,  reasonable  variations  from  this  profile  change 
the  thermal-noise  power  density  at  the  receiver  by  an  insignificant 
amount . 

Extensive  data  for  atmospheric  noise  at  the  earth's  surface  are 
available  fo'  various  seasons  and  geographic  locations.  VJe  use  data 


*As  shown  in  Sec.  IV,  the  transmission  properties  of  the  litho- 
sphere  depend  strongly  on  whether  abrupt  transitions  in  the  conductivity 
profile  are  present.  In  fact,  some  of  the  calculations  given  below  are 
for  the  Housley  BR  profile,  modified  to  have  a  sharp  transition  at  a 
depth  of  25  km. 
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compiled  by  Maxwell  (IOC?).  His  rms  noise-density  spectra  for  summer 
and  four  representative  locations  are  shown  in  Fig.  4.  Maxwell's  data 
indicate  that  these  rms  atmospheric  noise  values  are  exceeded  lO-to-20 
percent  of  the  time.  Thus,  a  system  designed  on  the  basis  of  the  rms 
atmospheric  noise  will  achieve,  or  exceed,  performance  specifications 
80-to-90  percent  of  the  time;  i.e.,  it  will  have  an  80-to-90  percent 
time  availability.  If  99-percent  time  availability  is  required,  then 
the  system  should  be  designed  on  the  basis  of  atmospheric-noise  spectral 
densities  about  10  dB  larger  than  those  shown  in  Fig.  4.  At  the  earth's 
surface,  the  electric  fields  associated  with  atmospheric  noise  in  the 
ELF/VLF/LF  bands  are  nearly  vertical;  whereas  the  magnetic  fields  are 
nearly  horizontal.  The  noise-field  components  at  depth  must  be  computed 
fo’'  each  conductivity  profile  from  the  surface-noise  data.  The  details 
of  this  calculation  are  given  in  Sec.  III. 


-120 
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III.  MATHEMATICAL  FORMULATION 

In  view  of  the  inherent  uncertainty  in  the  conductivity  depth- 
profiles  given  in  Sec.  II,  the  choice  of  computational  method  requires 
discussion.  It  is  tempting  to  argue  that  highly  simplified  analytic 
techniques  will  suffice,  since — regardless  of  the  mathematical  approach 
the  final  results  will  be  no  more  accurate  than  the  models  used  for  inputs. 
Although  clearly  correct,  this  argument  should  not  be  carried  too  far, 
lest  Important  general  characteristics  of  the  electromagnetic  fields  be 
obscured.  In  particular,  account  must  be  taken  of  the  depth  gradients 
of  conductivity  to  gain  insight  into  the  depth  dependence  of  the  signal 
and  atmospheric  noise  fields.  This  insight  is  needed  to  determine  the 
dependence  of  system  performance  on  antenna  burial  depths.  Clearly,  if 
either  transmitting  or  receiving  antennas  could  be  placed  at  relatively 
shallow  depths  above  the  waveguide,  system  costs  would  be  reduced  signifi¬ 
cantly.  Conductivity  gradients  also  strongly  affect  the  lateral  attenua¬ 
tion  in  the  waveguide.  Since  the  wavelengths  Involved  are  not  small 
compared  with  vertical  distances  over  which  the  electrical  properties  of 
the  lithosphere  change  substantially,  elkonal  methods  are  not  applicable 
and  full-wave  calculations  are  needed.  To  conipute  the  lateral  attenuation 
and  depth  dependence  of  the  fields,  we  have  used  a  method  that  accounts 
in  detail  for  the  vertical  inhomogeniety  of  the  lithospheric  conductivity. 
We  use  simpler  computational  methods  to  account  for  lateral  inhomogenieties 
and  to  determine  certain  system  performance  parameters. 

MODAL  SOLUTIONS 

We  wish  to  determine  the  spatial  dependence  of  the  electromagnetic 
fieldo  propagating  in  a  waveguide  that  is  strongly  inhomogeneous  in  the 
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vertical  direction.  The  curvature  of  the  earth  will  be  neglected  and 
horizontal  stratification  will  be  assumed  (see  Fig.  5).  The  method 
of  calculation  Is  a  version  of  the  waveguide  mode  analysis  described 
by  Budden  (1962).  This  method,  often  used  in  problems  involving  propa¬ 
gation  in  the  earth- ionosphere  waveguide  (e.g.j  Fields  1970),  has  been 
modified  in  this  study  for  applical ion  to  the  lithosphere. 

Although  the  best  antenna  configuration  for  use  in  a  lithospheric 
coninunication  channel  is  n't  known,  a  vertically  polarized  E-fleld 
antenna  is  a  likely  candidate.  We  thus  consider  plane  TM  wavegnlde 
modes  propagating  in  the  x-directlon,  for  which  the  electric  field,  E, 
and  the  magnetic  field  can  be  written 

i(tL)t-kn  S  x) 

E(x..)  .  ^  ' 


l(a)t“kn^S^x) 

V(x,z)  -  CyVyCz)  e 

In  Eq.  (1),  =  where  H  is  the  usual  magnetic  Intensity,  and 

f  are  the  electric  and  magnetic  permittivities  of  free  space,  w  -  Z^rf 
(where  f  is  the  wave  frequency),  t  is  tliie,  and  k  =  w/c  (where  c  is  the 
vacuum  speed  of  light).  MKS  units  will  be  used.  The  refractive  index, 
n,  is  given  by 


n^(z) 


e(z)  - 


iq(z) 


c 


(2) 


where  a  is  the  depth-dependent  conductivity  and  t  is  the  relative  electric 
permittivity.  The  quantity  n^  in  Eq.  (1)  is  the  refractive  index  at  the 
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depth,  z^,  of  minimum  conductivity  in  the  lithosphere.  can  be 

interpreted  as  the  sine  of  the  incidence  angle  at  the  depth  z^.  Incidence 
angles  at  depths  other  than  z^  can  be  determined  from  the  relation 

S(z)n(z)  “  n^S^  -  constant,  (3) 

which  follows  from  the  assumed  horizontal  stratification  and  is  essentially 
Snell's  law. 

The  relevant  Maxwell  equations  are: 

dE 

+  Ikn  S  1  -  -ikV  ; 

dz  o  o  J  y 


_JIL 

dz 


-ikn^(z)E^  ; 
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0  0 
n^(z) 


(6) 


Rather  than  solving  Eqs.  (A)  through  (6)  directly,  we  use  the  wave  ad 
mlttance.  A,  defined  by 


A  i  . 

and  the  related  quantity 


A- 
A+1  * 


(7) 


(8) 


By  combining  Eqs.  (4)  through  (8),  it  follows  that 
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dz 


Ik 
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n^(z)(W-l)^ 


n^(z) 


(W+1)^ 


(9) 


It  is  clear  from  Eq.  (1)  that  the  lateral  attenuation  in  the  waveguide 
Is  governed  by  the  product  n^S^.  Since  n^  Is  specified  once  a  lithospheric 
model  and  a  wave  frequency  have  been  selected,  only  must  be  calculated 
to  determine  the  attenuation.  This  calculation  Involves  the  solution  of 
Eq.  (9)  subject  to  the  appropriate  boundary  conditions,  whence  eigen¬ 
values  for  S  are  found.  A  guided  wave  In  the  lithosphere  can  produce 
o 

only  upgolng  waves  In  the  medium  above  the  earth's  surface,  since  this 
medium  Is  assumed  uniform  with  constant  Index  of  refraction,  n^^.  Thus, 
Immediately  above  the  earth's  surface,  the  wave  admittance  Is  given  by 


ACO"^)  - 


n. 


(10) 


For  propagation  under  the  ocean,  n^  Is  the  refractive  Index  of  sea  water, 
and  may  be  assumed  Infinite  for  the  frequencies  to  be  considered.  For 
this  case,  ACO"*")  ->■  «  and,  from  Eq.  (8),  the  boundary  condition  for  W  Is 

simply 


W(S  „,z-0)  «  1  , 

ON 


(lla) 


where  S  „  denotes  the  eigenvalue  of  S  for  the  mode.  For  the  case 
where  the  medium  above  the  earth's  surface  Is  free  space,  n^^  -  1;  and, 
using  Eqs.  (3),  (8),  and  (10),  the  boundary  condition  on  W  becomes 


W(S^j^,z.O) 


1/2 


1/2  ■ 


(11b) 
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n,e  conductivity  ptoflle.  glvnn  In  Snc.  II.  nnd  values  fat  t  given 
below,  ate  aulflelent  to  detemlne  n^t)  lot  the  various  models.  Once 
n^  1.  specified.  E,a.  (9)  and  (11)  comprise  a  closed  set  for  W  and 
These  coupled  equations  are  solved  readily  hy  straightforward  Iteration. 
Each  Iteration  requires  the  numerical  solution  of  Eq.  (9).  This  solution 
is  started  at  a  great  depth,  where  a  purely  downgoing  wave  Is  assumed 
as  an  Initial  condition.  In  each  solution,  several  starting  depths  and 
integration  step  sites  are  tried;  and  the  process  Is  truncated  when 
successive  trials  fall  within  a  prescribed  tolerance.  The  Integration 
is  terminated  at  the  earth's  surface.  For  a  profile  such  as  that  shown^ 
In  Fig.  2.  It  is  inconvenient  numerically  to  use  downgoing  waves  at 
great  depths  as  an  Initial  condition,  since  Integration  across  any  rones 
of  rapid  conductivity  change  (which,  for  the  case  shown  In  Fig.  2.  would 
be  at  14  km)  would  be  required  on  each  Iteration.  For  the  ptoflle  In 
Fig.  2.  we  found  that  the  use  of  W(r-14)  -  1  as  an  Initial  condition 
provided  good  accuracy  and  considerable  simplification.  This  Initial 
condition  is  equivalent  to  truncating  the  conductivity  profile  In  . 
perfect  conductor  at  z*-lA  km. 

Once  the  eigenvalue  Is  found  and  W  (and  hence  A)  Is  computed, 

the  attenuation  rate  for  the  h'"  mode  and  the  depth  dependence  of  the 
electromagnetic  fields  ate  easily  obtained.  The  attenuation  rate  Is 

given  by 


-  8.7  X  dB/km  ; 


(12) 


and,  by  Integrating  Maxwell's  equations,  it  follows  that 


Vy(z) 


exp 


[-  f  ■ 


(13) 


where  V  has  been  taken  to  equal  unity  at  the  "center"  of  the  waveguide; 

y 

l.e.  at  z“z  .  The  time  and  x-dependence  of  the  fields  have  been  sup- 
pressed  in  Eqs.  (13)  through  (15). 

ATMOSPHERIC  AND  THERMAL  NOISE 

To  calculate  achievable  transmission  ranges  and  data  rates  in  the 
lithosphere,  it  is  necessary  to  know  the  atmospheric  and  thermal-noise 
spectral  density  at  the  receiver  depth.  The  atmospheric  noise  fields 
at  depth  must  be  calculated  from  the  available  data,  which  are  typically 
in  the  form  of  vertical  electric  fields  measured  at  the  earth's  surface 
(see  Fig.  A).  The  quantities  that  relate  the  noise  fields  at  depth  to 
the  surface  noise  fields  will  be  called  noise-transfer  functions.  For 
uniform  media,  such  as  the  ocean,  the  transfer  functions  can  be  calculated 
easily  in  terms  of  plane  waves  and  the  well-known  Fresnel  reflection  and 
transmission  coefficients.  However,  when  the  conductivity  exhibits 
significant  depth-gradients  (e.g.,  as  in  Figs.  1-3),  the  transfer  functions 
cannot  be  computed  in  closed  form.  Gallawa  and  Paidle  (1972)  accounted 
for  the  vertical  Inhomogenlety  of  conductivity  by  representing  the  litho¬ 
sphere  by  several  uniform  layers  and  computing  the  transfer  function  by 
summing  the  absorptions  suffered  in  each  layer.  However,  even  this 
approach  is  not  adequate  for  the  problem  at  hand.  It  neglects  both  the 
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conversion  lose  suffered  at  the  earth's  surface  and  the  fact  that  the 
admittance  "seen"  by  a  dounuard  propagating  signal  Involves  an  Integral 
over  all  depths.  Also,  the  total  nolae  field  at  depth  should  not  he 
used,  instead,  the  component  of  the  noise  field  In  the  direction  of  the 
antenna  polarisation  Is  needed.  Although  the  noise  fields  above  the 
earth's  surface  are  nearly  vertically  polarized,  the  fields  In  the 
highly  refractive  earth  can  be  nearly  horizontally  polarized.  And  a 
vertically-oriented,  burled  electric  antenna  (a  likely  configuration) 

„111  be  1 . nsltlve  to  atmospheric  noise  than  would  a  horizontally 

oriented  one.  In  fact,  the  strength  and  polarization  of  the  subterranean 
nolae  fields  are  related  In  a  complicated  fashion  to  the  conductivity 
profile  and  change  continuously  with  depth.  A  full-wave  numerical 
calculation  of  these  fields  Is  thus  needed.  Fortunately,  the  procedures 
developed  to  solve  E,s.  (9)  and  (13)  through  (15)  can  be  used  essentially 

without  change. 

The  problem  at  hand  Is  to  relate  the  noise  fields  at  depth  to  the 
vertical  atmospheric  electric  noise  field.  E^.(0).  which  Is  assumed  known 
at  the  earth's  surface.  Atmospheric  noise  typically  propagates  above  the 
earth's  surface  as  a  TM  mode  with  an  Incidence  angle  larger  than.  say. 
70“.  Let  the  sine  of  this  angle  be  denoted  by  S^.  The  magnetic  noise 
field  at  the  surface  is  then  given  by 


(0)  -  • 


(16) 


The  depth-d.p«.d.nt  admittance  seen  by  this  signal  can  be  computed  by 
solving  Egs.  (8)  and  (9).  with  z^  -  0.  whence  and  n_,  -  1.  Note 

that  no  Iteration  1.  needed  since  8.  Is  specified,  and  only  a  single 


integration  of  Eq.  (9)  is  needed  for  each  profile  and  wave  frequency. 
Once  the  admittance  has  beer  thus  determined,  the  transfer  functions, 
Y^,  are  easily  computed  from  Eqs.  (13)  through  (15).  The  resulting 


equations  are: 


V  (z)  -  Y-E,^(0), 
"ya  1  za 


(17a) 


(17b) 


where 


(17c) 


1  S 


(18a) 


^2  -  '  2  "^1’ 
^  n^(z) 


(18b) 


Yj  -  Yj^/A(z). 


(18c) 


Strictly  speaking,  the  noise-transfer  functions  given  above  cannot  be 
completely  specified,  since  the  incidence  angle  (and,  hence,  S^) ,  is 
never  known  precisely.  However,  these  functions  are  very  insensitive 
to  reasonable  variations  in  S^.  For  eirample,  calculations  for  incidence 
angles  of  85°  and  70°  (a  realistic  range  of  values)  yielded  values  for 
Y^  that  differed  by  less  than  1  or  2  dB-an  insignificant  variation. 

The  results  given  in  Sec.  IV  are  computed  for  an  assumed  incidence  angle 
of  72°  (S  -  0.95),  and  should  very  accurately  approximate  the  transfer 
functions  for  all  obliquely  Incident  noise  fields. 


-25- 


For  propagation  above  the  earth's  surface  and  frequencies  in  the 
LF  band  and  below,  atmospheric  noise  is  always  much  stronger  than  thermal 
noise;  and  above-ground  long-wave  systems  are  atmospherically  noise 
limited.  However,  atmospheric  noise  is  attenuated  in  the  earth's  crust 
and  thus  decreases  w^th  increasing  depth;  and  thermal  noise  becomes 
greater  with  depth  because  of  the  higher  temperatures.  For  certain 
situations,  thermal  noise  could  be  the  limiting  factor  for  a  deeply 
buried  receiver.  The  thermal  noise  power  density,  N^,  will  be  computed 
from  the  simple  relation 

•=  1.38  X  lO'^^T  watts/Hz  ,  (19) 

where  T  is  the  lithospheric  temperature  at  the  receiver  depth  in  degrees 
Kelvin.  T  will  be  taken  from  the  Keller  model  as  shown  in  Fig.  2. 

Equation  (19)  for  is,  of  course,  oversimplified.  In  a  non-uniform 
medi^im  such  as  the  lithosphere,  the  antenna  receiver  "sees"  a  continuum 
of  temperatures  from  the  surrounding  material.  However,  the  error  in¬ 
curred  by  using  Eq.  (19)  is  quite  small,  as  can  be  seen  from  the  fact 
that  the  temperature  varies  by  less  than  a  factor  of  2  (less  than  3  dB) 
over  the  entire  range  of  depths  shown  in  Fig.  2.  Further,  other  avail¬ 
able  lithospheric  models  show  temperatures  that  differ  from  those  shown 
in  Fig.  2  by  no  nwre  than  20  percent.  Gallaua  and  Haidle  (19^  <i)  have 
carried  out  detailed  calculations  of  N.^.,  taking  proper  account  of  antenna 
directivity  and  Integrating  over  the  contribution  from  the  surrounding 
medium.  As  would  be  expected  from  the  above  discussion,  their  results 
differ  from  those  obtained  from  Eq.  (19)  by  a  few  dB,  at  most.  A  few 
dB  are  insignificant  compared  with  other  factors  that  influence  expected 


system  performance. 
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rate  of  R  bits  per  second  is  given  by 


(20) 


req 


E,  IS  the  enets,  pet  hit  et  the  teeeivet.  end  ie  the  etteettve 

(l.e..  the  post-preeeeelng)  nolee-pewet  density.  We  assume 
vhleh  is  sufficient  to  achieve  an  .trot  rate  of  10  with  simp 
lation.  Equation  (20)  thus  becomes 


S  =  ^  ’ 

req  etl 


(21) 


Pot  ahove-gtound  eo^unlc.tlons .  It  Is  cuatomat,  to  Include  a  signal 
»tgln  to  insure  tellahlllty  even  during  adverse  fluctuations  In  propaga¬ 
tion  conditions.  Such  a  margin  has  not  been  Included  In  Eg.  (21),  because 
(unllhe  the  Ionosphere)  the  lithospheric  waveguide  should  he  a  very  stable 

communication  channel. 

in  the  following  discussion,  several  Idealizations  are  made  to  render 
the  calculations  of  the  received  signal  and  noise  trectahle.  These  approxi¬ 
mations  are  necessary  because  the  problem  of  calculating  the  directivity 
and  effective  area  of  an  antenna  Imbedded  In  a  medium  that  Is  both  con- 
ductive  and  inhomogeneous  has  never  been  solved.  The  directivity 
of  the  receiving  antenna  Is  needed  because  the  signal  end  atmospheric  noise 
utrlv.  from  very  different  dlrectlons-the  signal  propagation  being  more 
ut  less  lateral,  whereas  the  atmospheric  noise  propagates  mainly  vertically 
in  the  earth's  crust.  «e  proceed  with  the  c.lculatlon  as  If  the  receiv¬ 
ing  antenna  were  Imbedded  In  a  uniform  medium.  Although  difficult  to 


justify  rigorously,  this  approach  should  give  a  reasonably  good  estimate 
of  the  required  power.  For  all  of  the  lithospheric  models  except  one 
(Levin  "wet"  model.  Fig.  1),  the  refractive  index  in  the  crust  is  fairly 
uniform  and  has  a  small  imaginary  part  (implying  a  dielectric-like 
medium)  at  depths  where  the  conductivity  is  small  and,  hence,  where  a 
receiver  is  apt  to  be  placed.  Moreover,  although  inaccurate  in  detail, 
the  expressions  used  do  incorporate  the  Important  features  of  the 
propagation  and  reception.  For  example,  the  expression  used  for  antenna 
directivity  does  correctly  account  for  the  fact  that  a  vertical  E-fleld 
antenna  is  sensitive  only  to  vertical  electric  fields,  and  thus  mitigates 
much  of  the  atmospheric  noise  that  is  mainly  horizontally  polarized  at 
depth.  The  results  thus  should  be  adequate  for  this  feasibility  study. 

For  illustrative  puvposes,  we  assume  that  the  receiving  antenna  is 
a  short,  vertical  electric  dipole.  This  configuration  is  a  logical 
choice,  since  it  should  fit  easily  into  a  narrow,  vertical  borehole, 
and  will  be  relatively  insensitive  to  downward  propagating  atmospheric 
noise.  Note  that  detailed  antenna  design  is  beyond  the  scope  of  this 
study,  and  future  study  could  well  demonstrate  that  some  other  configura¬ 
tion  is  more  attractive.  However,  it  is  unlikely  that  alterations  in 
receiver  design  would  change  the  forthcoming  conclusions  regarding  system 
feasibility.  For  example,  any  realistic  antenna  must  be  fairly  small 
electrically  because  of  the  large  wavelengths  in  the  crust  (>  1-lOkm) , 
and  all  small  antennas  have  about  the  same  gain.  Further,  as  is  evident 
from  the  numerical  results  given  in  Sec.  IV,  system  feasibility  is 
dominated  by  the  attenuation  of  the  signal  and  atmospheric  noise  in  the 
crust.  Neither  of  these  factors  depend  in  any  way  on  hardware  design. 


For  the  effective  area  of  the  receiving  antenna,  we  use 


3/2  IT  sin  9 
2 

0)  y  e  e  (z) 
o  o 


(meters)' 


(22) 


where  the  angle,  6,  is  measured  relative  to  the  vertical.  The  incident 

intensities  I  and  I  ,  associated  with  the  signal  and  atmospheric  noise 
s  a 

density,  respectively,  are  given  by 


and 


2  2 

I  -  1/2  VTTy  Re(n(z))  E  (z)  watts/m  , 

s  ’  o  o 

I  -  1/2  Ve  /p  Re(n(z))  E^(z)  watts/m^-Hz  , 

a  •  O  a 


(23) 


(2A) 


where  E  is  the  total  signal  electric  field  and  E^  is  the  noise-field 
spectral  density.  Noting  that  E^  sin^0  =  E^,  etc.,  it  follows  that 

Re(n(z))  E^(z) 

Signal  Power  =  - 2 -  watts  ,  (25) 

160  k  e(z) 

and 

,  Re(n(z))  E^„(z) 

Atmospheric  Noise _ za - watts/Hz  ;  (26) 

Power  Density  "  ^2^^^^ 

whereas  the  internal  noise  density  is  simply 

N^  -  1.38  X  lO"^^  TF  watts/Hz  ,  (27) 

where  F  is  the  receiver-noise  figure.  The  effective  noise  density  at 
the  receiver  is  thus 
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Re(n(z))  (z) 

N  -  - - - — -  +  1.38  X  10  TF  watts/Hz  ,  (28) 

160  k^e(z)G 

P 


where  G  is  the  processing  gain;  i.e.,  the  factor  by  which  the  effective 
P 

atmospheric  noise  can  be  suppressed  by  using  non-linear  elements  in  the 
receiver  circuit. 

By  Inserting  Eqs.  (25)  and  (28)  into  Eq.  (21),  we  find  the  following 
expression  for  the  required  vertical  electric  field  at  the  receiver 


(E^z)) 


req 


^za^^^  (160)  (1.38  X  10~^^)k\  (z) 

G  Re(n(z)) 

P 


(29) 


I  The  receiver  Is  assumed  to  be  at  a  depth,  z.  The  final  step  for  computing 

required  power  Is  to  relate  the  signal  electric  field,  E^(z),  to  the 

i 

'  transmitter  power.  This  can  be  done  approximately  by  noting  that  the 

I  waveguide  signal  spreads  nearly  cyllndrlcally  from  the  source  and  also 

suffers  an  exponential  attenuation  given  by  F.q.  (1?).  At  a  distance,  d, 
from  the  receiver,  the  signal-power  density  can  be  written  In  the  form 

i 

f 

nP  exp(-10  d/4.3]  „ 

I  - - - -  ■  ■  - -  watts/m  .  (30) 

s  Zrrd  h 

eff 

Only  a  single  waveguide  mode  (the  N*"^)  has  been  used  In  Eq.  (30),  since 
we  are  considering  distances  large  enough  that  only  the  least  attenuated 
mode  will  contribute  significantly.  is  the  total  transmitter  power, 
and  n  is  the  transmitter  efficiency;  viz.,  the  fraction  of  the  total 
power  that  Is  radiated  Into  the  dominant  (N*"^)  waveguide  mode.  The 
quantity  Is  defined  by 


I 


,0  ^  Re(n(z))  E^(.) 

•'-<c  Re(n(z^))  E^(z^) 

and  represents  the  nominal  width  of  the  waveguide  region  that  contains 
.„st  o£  th«  .Ignal  povar.  From  E,s.  (17b).  (23).  and  (29)  through  (31). 
it  follova  that  £or  a  racalvar  locatad  at  the  "uenter”  (t-r„)  of  the  wave- 
guide,  the  required  transmitter  power  is  given  by 


67iRh^,,d  fRe(n^)Y2(z^)  eJ^(O) 

(Po^  "  “i; - 1207^ 

req  P 


/4k^c  ( 

[— 


Z  )  \  _0'i 

-2-  1.38x10  T(z  )F 


[10'\j^d 


receiver  depths  other  than  z^,  the  required  power  can  be  estimated 


from: 


(P(t))_-D(z)(P„)  , 


where  the  depth  degradation  factor,  D,  is  given  by 


Re(n(z))Y?(z)E2  (0) 


/  *2 

E^(z„)  I  120’  g, 

eIm  'Re(n(t„))Y^(a^)E^,(0) 
^  120ir  G 

P 


1.38xlO“^^T(z)F 


38xlO“^^T(z^)F 


where  z  denotes  the  receiver  depth,  the  transmitter  depth  being  unspecified. 
For  systems  limited  by  atmospheric  noise,  D  simplifies  to 
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D(z) 


E^(z  )  Re(n(z))  Y^Cz) 

A,  Z  O _ __± - 

^  E^(z)  Re(n(z^))  Y^(z^) 


(35) 


For  systems  limited  by  thermal  noise: 


E^(z  >  eiz)  T(z) 

D(z)  %  ^ -  • 

e2(z)  e(z^)  T(7^) 

The  rationale  used  in  selecting  the  various  input  parameters 
in  Sec.  IV. 


(36) 

is  discussed 


! 
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IV.  NUMERICAL  RESULTS  AND  DISCUSSION 

Pacific-Sierra's  long-wave  propagation  code  is  used  to  calculate 
the  attenuation  rates,  and  the  signal-  and  noise-wave  functions,  for 
several  of  the  model  waveguides  shown  in  Figs.  1  through  3.  The  results 
of  these  calculations  are  given  in  this  section.  Detailed  results  are 
also  given  for  the  required  power  as  a  function  of  data  rate,  trans 
mission  range,  and  receiver  burial  depth.  In  all  cases,  the  relative 
electric  permittivity,  e,  is  assumed  to  have  a  constant  value  of  ten. 

ATTENUATION  RATES 

Figure  6  shows  the  calculated  attenuation  rates  corresponding  to 
the  various  model  conductivity  profiles.  Results  are  shown  only  for  the 
least-attenuated  waveguide  inodes;  viz.,  the  TEM  mode  for  'requencies 
lower  than  3  or  A  kHz,  and  the  lowest  order  TM  mode  for  higher  frequencies. 
The  results  showed  that  Includion  of  only  these  moda'a  is  adequate  for 
determining  maximum  transmission  ranges.  The  decision  to  restrict 
attention  to  frequencies  between  100  Hz  and  100  kHz  was  somewhat,  but 
not  entirely,  arbitrary.  It  is  very  difficult  to  radiate  useful 
amounts  of  power  at  frequencies  as  low  as  100  Hz.  At  frequencies 
higher  than  about  100  kHz,  the  computational  methods  used  here  become 
inconvenient  due  to  the  necessity  of  t.  -.aining  large  numbers  of  wave¬ 
guide  modes  rather  than  the  one  or  two  least-attenuated  ones.  Ray¬ 
tracing  methods  are  probably  preferable  to  modal  analysis  at  these 
higher  frequencies.  Also,  the  degrading  effects  of  scattering  from 
irregularities  in  the  crust  would  be  expected  to  become  more  pronounced. 
Nonetheless,  it  could  well  prove  worthwhile  to  extend  consideration 
to  frequencies  outside  of  the  ELF/VLF/LF  bands  treated  here. 


^‘*tenuation  Rate  vs.  Frequency  ^cr  Four  *^del  Conductivity  Profiles 
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The  results  shjwn  in  Fig.  6  apply  to  propagation  in  the  continental 
lithosphere,  since  a  transition  to  free  space  was  assumed  at  z^O  (Fig.  5) 
Attenuation  rates  were  also  calculated  for  the  case  where  the  crust 
interfaces  with  sea  water  (approximated  by  a  perfect  conductor)  at  z'O. 
These  results,  which  will  not  be  presented  in  detail,  are  virtually 
identical  to  those  s’.iown  in  Fig.  6  for  frequencies  higher  than  a  few 

it 

hundred  Hertz.  At  lower  frequencies,  the  effect  of  assuming  a  sea-water 
"top"  on  the  waveguide  la  to  reduce  the  attenuation  from  the  values  shown 
particularly  for  the  Keller  and  Houaely  BR  models.  This  behavior  is  to 
be  expected  since,  for  extremely  low  frequencies  where  the  skin  depth 
is  large,  sea  water  reflects  energy  that  would  otherwise  leak  out  of  the 
waveguide . 

The  general  dependence  of  attenuation  on  frequency,  as  shown  in 

Fig.  6,  is  quite  similar  to  that  observed  for  long-wave  propagation  in 

the  earth-ionosphere  waveguide;  viz.,  relatively  low  attenuation  at 

frequencies  less  than  a  few  hundred  Hertz  or  higher  than  10  kHz,  and  a 

"forbidden"  band  in  the  l-to-5  kHz  portion  of  the  spectrum.  However,  in 

the  earth-ionosphere  waveguide,  the  attenuation  rate  is  typically  2x10 
_3 

to  4x10  dB/km  in  the  VLF  band.  And  the  attenuation  seldom  exceeds 
_2 

10  dB/km,  even  under  disturbed  ionospheric  conditions.  By  comparison, 
the  computed  attenuations  in  the  lithospheric  waveguide  are  an  order  of 
magnitude  or  more  larger. 

The  fact  that  the  attenuation  depends  strongly  on  the  minimum 
value  of  conductivity  and  the  conductivity  depth-gradients  is  clearly 
illustrated  by  the  results  shown  in  Fig.  6.  The  minimum  conductivity  is 

* 

The  reader  is  cautioned  against  applying  these  results  to  propaga¬ 
tion  under  the  oceans,  however,  since  the  models  shown  in  Figs.  1  through 
3  were  developed  to  represent  the  continental  lithosphere. 
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about  10  ^mhos/m  for  the  L.»vln  (wet)  model,  about  10  ^mhos/m  for  both 

-9 

the  Levin  (dry)  and  Keller  models,  and  about  5x10  mhos/m  for  the  Housley 
BR  model.  As  would  be  expected,  the  Housley  BR  model  yielded  the  lowest 
attenuation,  whereas  that  yielded  by  the  Levin  (wet)  model  was  by  far  the 
highest.  However,  the  attenuation  rates  for  the  Levin  (dry)  model  and 
the  Keller  model  are  quite  different,  even  though  the  minimum  conductivity 
Is  nearly  the  same  for  the  two  models-  At  a  frequency  of  1  kHz,  the 
calculated  attenuation  rate  Is  about  0.18  dB/km  for  the  Keller  model  and 
about  0.28  dB/km  for  the  Levin  (dry)  model.  The  more  favorable  attenuation 
In  the  Keller  model  of  the  lithospheric  waveguide  Is  due  In  part  to  the 
assumed  presence  of  a  rapid  transition  In  conductivity  at  a  depth  of 
about  14  km  (Fig.  2) . 

It  Is  Interesting  to  ccupare  the  results  shown  In  Fig.  6  with  the 
attenuation  rates  computed  by  Spie"^  and  Wait  (1971)  for  a  step-function 
model  of  the  conductivity  depth-profile.  For  an  assumed  waveguide 
conductivity  of  10  ^  mhos/m,  their  calculated  attenuation  rates  are 
typically  a  factor  of  two  or  three  lower  than  those  shown  in  Fig.  6 
for  the  Keller  and/or  Levin  (dry)  models.  Thi»se  lower  attenuation  rates 
are  due  largely  to  Spies  and  Walt's  use  of  sharp,  highly  reflecting  upper 
and  lover  boundaries  rather  than  the  gradual  "boundaries"  used  here.  Of 
course,  the  conductivities  used  by  Spies  and  Walt  In  their  step-function 
representation  of  the  waveguide  are  Intended  to  be  "effective"  values, 
averaged  over  some  appropriate  range  of  depths,  rather  than  minimum  values. 
In  this  context,  the  effective  conductivity  of  the  Keller-model  waveguide 
would  be,  say,  5  or  6x10  ^mhos/m,  even  though  the  minimum  conductivity 
is  about  10  ^mhos/m. 
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The  significant  effect  of  conductivity  depth-gradients,  and  the 
presence  or  absence  of  sharp  transitions  in  the  conductivity  depth- 
profile,  can  be  illustrated  by  examining  the  llousley  BR  model  more 
closely.  For  reasons  discussed  in  Secs.  II  and  III,  the  results 
shown  in  Fig.  6  are  not  computed  for  the  Housley  BR  model  precisely 
as  shown  in  Fig.  3  (p.  13).  Instead,  the  conductivity  values  shown 
are  used  for  depths  down  to  25  km,  where  a  sharp  transition  to  a  highly 
conducting  sub-basement  (i.e.,  a  highly  reflecting  "boundary")  is 
assumed.  As  a  check  on  sensitivity,  additional  calculations  of  the 
attenuation  rates  for  the  Housley  BR  model  are  carried  out  for  assumed 
lower  "boundary"  depths  other  than  25  km.  Sample  results  are  given 
in  Fig.  7  for  a  frequency  of  1  kHz.  The  attenuation  is  seen  to 
increase  markedly  as  the  assumed  layer  depth  is  increased.  Results 
corresponding  to  the  Housley  BR  profile,  precisely  as  shown  in  Fig.  3, 
can  be  obtained  by  taking  the  depth  of  the  assumed  "boundary"  to  be 
infinite.  For  this  case  (see  Fig.  7),  the  attenuation  rate  is  about 

it 

twice  as  large  as  when  the  profile  is  terminated  at  25  km,  and  is 
nearly  as  large  as  for  the  Keller  model  (see  Fig.  6).  Of  course, 
if  the  assumed  transition  at  14-km  depth  were  removed  from  the  Keller 
model,  or  if  it  were  assumed  to  occur  at  a  greater  depth,  then  the 
calculated  attenuation  rates  would  be  larger  than  shown  in  Fig.  6. 


*It  should  not  be  inferred  from  Fig.  7  that  the  plac^ent  of  the 
reflecting  layer  at  a  depth  of  25  km  causes  the  results  shown  in  Fig.  6 
(and  to  be  shown  below  in  figures  for  the  Housley  BR  model)  to  be  unduly 
optimistic.  It  can  be  inferred  from  seismic  data  that,  for  many  regions, 
a  highly  reflecting  transition  exists  at  depths  shallower  than  25  km,  in 
which  case  the  Housley  BR  results  would  be  more  favorable  than  those  on 
whicn  the  conclusions  of  this  report  are  based. 
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SIGNAL  AND  NOISE  FIELD-STRENGTH  PROFILES 

Equations  (13)  through  (15)  have  been  integrated  to  obtain  the 
field-strength  profiles  associated  with  the  least  attenuated  waveguide 
modes.  Figures  8  and  9  show  sample  results,  illustrating  the  calculated 
field-strength  profiles  for  the  Keller  model  and  frequencies  of  1  kHz 
(TEM  mode)  and  10  kHz  (TM  mode).  Recall  that  in  all  cases  the  fields 
are  normalized  such  that  V-  1  at  the  "center"  of  the  waveguide;  l.e.. 
at  the  depth  of  minimum  conductivity.  Thus,  for  the  Keller  model,  '*( 
has  been  set  equal  to  unity  at  a  depth  of  about  10  km. 

One  of  the  goals  of  this  study  is  to  estimate  the  dependence  of 

expected  system  performance  on  the  burial  depth  of  the  receiving  antenna. 

Thus,  for  a  vertically  oriented  E-field  receiver— the  configuration 

assumed  for  illustrative  purposes— the  depth  dependence  of  IpJ  is  of 

particular  interest.  Figures  8  and  9  show  that  |eJ  has  a  broad 

maximum  centered  at  a  depth  of  about  10  km.  However,  for  depths  shallower 

than  5  or  6  km,  Ie  I  decreases  rapidly  as  the  depth  is  decreased.  This 
z  ' 

behavior,  which  is  characteristic  of  a  trapped  mode  in  a  waveguide 
"centered"  at  a  depth  of  10  km,  is  due  to  four  physical  effects: 

1)  energy  is  absorbed  as  the  signal  propagates  upwards  toward  the 
earth's  surface;  2)  the  conductivity  and,  hence,  the  admittance  of 
the  medium  increases  as  the  depth  is  reduced— causing  the  electric 
fields  to  be  reduced  relative  to  the  magnetic  field;  3)  the  signal 
suffers  gradient  reflection  as  it  "tries"  to  propagate  upwards;  4) 
since  the  refractive  index  Increases  as  the  depth  is  reduced,  the 
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"wave"  normal  becomes  more  nearly  vertical  as  the  signal  propagates 
upward.*  This  effect  causes  |E^j  to  be  reduced  relative  to  both  |E^1 
and  I/*/"  |.  For  depths  greater  than  about  4  km  and  the  two  frequencies 
corresponding  to  Figs.  8  and  9,  greater  than  which 

indicates  a  signal  traveling  mainly  in  a  horizontal  direction.  At 
depths  shallower  than  4  km,  the  propagation  direction  has  become  more 
nearly  verticil  than  horizontal.  Of  the  four  effects  described  above, 
absorption  becomes  relatively  more  important  as  the  frequency  is 
Increased,  since  the  skin  depths  are  reduced,  as  are  the  admittances 
and  refractive  indices. 

Figure  10  shows  the  frequency  dependence  of  the  relative  vertical 
electric  field  strength  o  several  depths  for  the  Keller  model.  As 
indicated  on  the  figure,  the  relative  field  strength  is  normalized  to 
unity  at  a  dep .  of  10  km.  As  would  be  expected,  the  relative  field 
strength  becomes  smaller  as  the  depth  becomes  shallower.  Field-strength 
profiles  have  also  been  calculated  for  the  Housely  and  Levin  models, 
but  in  the  Interest  of  brevity  are  not  presented  in  detail.  For  depths 
shallower  than  about  10  km,  the  field-strength  profiles  based  on  the 
Housley  BR  model  are  very  nearly  the  same  as  those  shown  for  the  Keller 
profile.  This  similarity  is  to  be  expected,  of  course,  since  the 
conductivity  profiles  in  the  two  models  are  identical  for  most  of  the 
O-to-10  km  depth  range.  On  the  other  hand,  the  Levin  (dry)  model  is 

*This  interpretation  is  intended  solely  to  provide  an  intuitive 
understanding  of  the  numerical  results  and  should  not  be  taken  too 
literally.  Strictly  speaking,  of  course,  the  concepts  of  wave  normal 
and  propagation  direction  do  not  apply  when  the  medium  has  large  spatial 
gradients,  as  is  the  case  here. 
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seen  (Fig.  1,  p.  9)  to  be  much  more  highly  conductive  than  either  the 
Keller  or  the  Housley  models  for  depths  shallower  than  about  8  km. 
Consequently,  for  depths  shallower  than  8  km,  the  relative  field 
strengths  for  the  Levin  (dry)  profile  are  orders  of  magnitude  smaller 
than  those  shown  In  Figs.  8  through  10. 

Noise-transfer  functions  have  been  computed  from  Eq.  (18). 

Figures  11  and  12  show  ]Y2|  and  |y^|  as  a  function  of  depth  for  the 
Keller  model,  and  frequencies  of  1  kHz  and  10  kHz.  Recall  that  the 
surface-noise  fields  shown  In  Fig.  4  should  be  multiplied  by  [y^I  or 
|y^|  to  obtain  the  vertical  or  horizontal  electrlc-nolse  fields, 
respectively,  at  a  given  depth.  The  noise  propagates  downward  from 
the  earth's  surface,  whereas  the  signal  propagates  upward  (and  laterally) 
from  the  center  of  the  waveguide.  Thus,  absorption  and  reflection  tend 
to  reduce  IY2I  and  |y^|  as  the  depth  Is  Increased.  This  behavior  Is 
different  than  noted  above  for  the  signal  fields,  which  suffer  reflec¬ 
tion  and  reflection  losses  as  the  depth  Is  decreased.  These  state¬ 
ments  do  not  Imply  that  the  transfer  functions  must  necessarily  decrease 
monotonlcally  with  Increasing  depth.  As  Illustrated  by  Fig.  11,  for 
example,  refractive  effects  can  be  dominant  at  the  lower  frequencies 
and  actually  cause  Iy^I  to  have  a  larger  value  near  the  center  of  the 
waveguide  than  near  the  earth's  surface.  Further,  since  the  conduc¬ 
tivity  becomes  low  near  the  center  of  the  waveguide,  the  crust  becomes 
quite  transparent  at  these  depths  and  standing  wave  patterns  In  the 
transfer  functions  can  occur.  An  example  of  such  a  pattern  Is  shown 
In  Fig.  13. 


Noise-lransfer  functions 
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Figures  11  and  12  show  that  [Y^!  is  considerably  smaller  than 
lY^l  for  most  depths  of  interest.  This  behavior  was  found  to  occur 
for  essentially  all  of  the  models  and  frequencies  for  which  calculations 
were  made.  Thus,  at  least  for  the  models  used  here,  the  vertical  noise 
electric-field  component  is  typically  mu^'h  smaller  than  the  horizontal 
one;  and  a  vertically  oriented  electric  antenna  should  therefore  be  less 
severely  degraded  by  atmospheric  noise  at  depth  than  a  horizontally 
oriented  one.  This  distinction  has  often  been  overlooked  (e.g.,  Gallava 
and  Haidle,  1972),  and  occurs  simply  because  low -frequency  electromag¬ 
netic  waves  tend  to  be  refracted  strongly  toward  the  vertical  upon 

entering  the  earth  from  free  space. 

The  frequency  dependence  of  the  noise-transfer  function  of  interest, 
Iy  1 .  is  shown  in  Fig.  13  for  the  Keller  model  at  several  depths.  At 

I  'I  I  * 

the  lowest  frequencies  considered,  the  vertical  component  of  atmospheric 
noise  tends  to  increase  with  increasing  depth,  indicating  a  dominance 
of  refraction  over  absorption.  Thus,  somewhat  surprisingly,  the  vertical 
component  of  atmospheric  noise  at  ELF  is  found  to  be  stronger  at  a  depth 
of  10  km  than  at  a  depth  of,  say,  2  km.  Of  course  the  total  normalized 
noise  power  decreases  more-or-less  monotonically  with  increasing  depth. 

At  the  higher  frequencies,  IY2I  is  seen  to  decrease  ar  the  depth  is 
increased,  indicating  that  absorption  is  the  dominant  mechanism 
determining  field  strength. 

DEPENDENCE  OF  SIGKAL-TO-NOISE  RATIO  ON  RECEIVER  DEPTH 

We  next  consider  the  depth-degradation  factor,  D(z)  ,  which  is 
simply  the  ratio  of  the  signal-to-noise  ratio  at  a  given  depth  to  that 


Noise- tronsfer  function,  IY2I 
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at  the  nominal  center  of  the  waveguide.  Thus,  In  the  Keller  model, 
for  example,  D(z)  Is  the  power  required  to  transmit  to  a  receiver  at 
a  depth,  z,  divided  by  the  power  required  to  transmit  to  a  similar 
receiver  burled  at  a  depth  of  10  km;  l.e.,  D(z)  Is  the  factor  by  which 
the  transmitter  power  must  be  changed  If  the  receiver  Is  moved  from 
a  10-km  depth  (Keller  model)  to  some  other  depth.  Since  10-km  deep 
holes  are  very  expensive  to  bore,  we  are  particularly  Interested  In 
determining  the  additional  power  that  would  be  required  If  the  receiver 
were  placed  at  some  lesser  depth. 

For  all  of  the  models  considered  In  this  report,  atmospheric  noise 

* 

dominates  the  thermal  noise  for  all  depths  and  frequencies  of  Interest. 
Thus,  the  simplified  Eq.  (35)  may  be  used,  and  D(z)  expressed  solely  as 
a  function  of  the  signal  field-strength  profiles  and  Some  examples 

of  the  depth-degradation  factor  are  given  for  the  Keller  model  (1  and 
10  kHz)  and  the  l.evln  (dry)  model  (1  kHz)  In  Fig.  14.  For  the  Keller 
model,  which  exhibits  a  fairly  constant  conductivity  depth-gradient 
between  0  and  10  km,  relatively  little  degradation  would  be  suffered 
by  reducing  the  receiver  depth  from  10  km  to,  say,  5  or  6  km.  Moreover, 
as  discussed  below,  the  penalty  In  transmission  range  Is  remarkably 
moderate  even  for  receiver  depths  as  shallow  as  3  km.  As  Indicated 
earlier,  the  Levin  (dry)  model  exhibits  a  rather  large  and  abrupt 

*Thls  dominance  occurs  because,  even  after  allowing  for  losses 
suffered  In  propagating  from  the  earth's  surface  to  the  nominal  center 
of  the  waveguide,  the  atmospherlc-nolse  density  was  found  to  be  stronger 
than  the  thermal-noise  density.  It  Is,  of  course,  possible  to  envision 
realistic  situations  In  which  the  atmospheric  noise  would  be  so  highly 
attenuated  as  to  be  no  longer  dominant.  A  case  In  point  would  be  If 
the  receiver  were  situated  In  the  suboceanlc  lithosphere. 


Degradation  Factor  vs.  Receiver  Depth 
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increase  In  conducrlvlty  as  Che  deprh  Is  reduced  ro  less  chan  abour 
8  k«.*  A  major  result  of  this  assumed  transition  In  conductivity  Is 
lllustra-.ed  in  Fig.  U:  D(z)  becomes  so  large  as  to  essentially  rule 
out  any  possibility  of  using  a  receiver  depth  shallower  than  7  or  8  hm. 
Thus,  to  the  extent  that  borehole  depth  affects  total  system  cost,  the 
attractiveness  of  a  lithospheric  communication  system  depends  quite 
strongly  on  the  gradient  of  conductivity  in  the  crust  at  depths  between 

0  and  10  km. 

Figure  15  shows  the  frequency  dependence  of  D(z)  for  the  Keller 
model  and  several  depths  of  Interest.  Of  course,  0(0-1  at  a  depth  of 
10  km.  AS  might  be  expected,  higher  frequency  signals  are  much  more 
sensitive  to  receiver  burial  depth  than  are  lower  fre  .uency  ones.  The 
large  Increase  in  0(0  for  shallow  depths  exhibited  at  frequencies 
higher  than  10  kHz  Is  due  to  two  factors:  1)  the  signal  gets  weaker 
as  the  surface  of  the  earth  Is  approached,  due  to  absorption  suffered 
in  propagating  upward  from  near  the  center  of  the  waveguide;  2)  the 
atmospheric  noise  gets  stronger  near  the  earth's  surface,  since  the 
absorption  ft  suffers  In  propagating  downward  Is  reduced.  For  reasons 
given  above.  D(z)  for  the  Housley  BR  model  Is  very  nearly  the  same  as 
that  shown  In  Fig.  15.  whereas  that  for  the  Levin  model  is  prohibitively 

large  at  depths  shallower  than  7  or  8  km. 

tdfai.;mt.;<;TQN  RANGES  AND  POWER  REQUIREMENTS 

The  dependence  of  required  power  on  data  rate,  transmission  range, 
transmitter  efficiency,  frequency,  and  receiver  burial  depth  has  been 

NhlTbehavlor  would  be  expected  In  regions  having  a  thick,  highly 
conducting,  sedimentary  layer. 
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calculated  from  Eqs.  32  and  33.  Values  for  rms  atmospheric  noise  corres 
ponding  CO  Colorado  summer  are  taken  from  Fig.  4,  (p.  15).  As  discussed 
in  Sec.  II,  the  use  of  rms  noise  spectra  implies  a  system  time  avail¬ 
ability  of  80-to-90  percent.  Also,  we  assume  no  processing  gain,  and 
have  thus  taken  Gp  equal  to  unity.  This  assumption  is  probably  quite 
conservative,  since  lO-to-20  dB  of  noise  suppression  can  typically  be 
achieved  in  above-ground  VLF/ELF  systems  {Evans  and  Griffiths^  1972). 

If,  say,  10  dB  of  noise  suppression  can  be  achltved  in  a  lithospheric 
system,  then  the  power  requirements  given  below  would  suffice  for  a 
time  availability  approaching  99  percent,  rather  than  the  above-stated 
80-to-90  percent.  For  purposes  of  presentation,  it  is  convenient  to 
define  a  normalized,  required  effective  power  (REP)  as  follows: 


n  P 

REP-— ^  (Watts), 


(37) 


where  P^^^  is  the  total  power  that  must  be  delivered  by  the  transmitter,* 
as  computed  from  Eq.  (33).  The  transmitter  efficiency,  n,  denotes  total 
efficiency;  viz.,  the  ratio  of  the  power  radiated  into  the  least- 


attenuated  mode  to  the  total  system  power;  n  thus  accounts  for  antenna 
radiation  efficiency,  the  excitation  efficiency  of  the  dominant  wave¬ 
guide  mode,  transmission  line  losses,  etc.  Therefore,  REP  denotes 


* 

Strictly  speaking,  REP  has  units  of  Joules/bit,  rather  than  Watts, 
since  R  has  units  of  blts-per-sec,  and  n  is  dimensionless.  The  distinc¬ 
tion  vanishes  for  a  data  rate  of  1  bps,  since  the  energy  per  bit  is  then 
numerically  equal  to  the  power.  A  more  precise  definition  than  given 
by  Eq.  37  would  be 


REP 


P  R 
_ o 

R 


P 

req 


where  Rg-l  bit-per-second.  In  the  following  discussion  and  figures,  the 
factor  Rq  is  suppressed,  and  REP  will  be  expressed  as  the  effective  power 
required  to  transmit  1  bit-per-second. 


the  power  that  must  be  radiated  Into  the  dominant  mode  to  achieve  a 
data  rate  of  1  bit-per-second  (bps).  REP  is  a  convenient  parameter, 
since  it  is  independent  of  n  (see  Eq.  32).  which  is  extremely  difficult 
to  calculate  for  a  transmitter  in  an  inhomogeneous  conducting  medium. 

It  is,  of  course,  a  simple  matter  to  compute  the  required  power, 
from  REP,  once  n  has  been  calculated  and  a  data  rate,  R,  chosen. 

The  transmitter  efficiency,  n,  cannot  now  be  specified  accurately. 
However,  some  assumption  must  be  made  regarding  achievable  transmitter 
efficiencies  if  system  feasibility  and  attractiveness  are  to  be  assessed 
quantitatively.  To  undertake  such  an  assessment,  we  assume,  without 
delving  into  the  specific  transmitter  configuration,  an  efficiency 
that,  based  upon  present  understanding,  does  not  appear  to  be  unduly 
optimistic.  Specifically,  in  addition  to  detailed  calculations  of  REP, 
results  will  be  given  for  P^^^  for  the  case  wl.ei  e  n/R  =  10  ;  i.e., 

where  the  total  required  power  exceeds  REP  by  a  factor  of  1000.  This 
choice  is  by  no  means  arbitrary.  Gallam  and  Haidle  (1972)  have 
estimated  that  efficiencies  of  a  few  percent  should  be  achievable  for 
a  VLF  transmitter  utilizing  a  deeply  buried  vertical  antenna.  Thus, 
at  VLF,  our  choice  of  n/R“10"^  can  be  interpreted  as  corresponding  to 
an  efficiency  of  a  few  percent  and  a  data  rate  of  a  few  tens  of  bps. 

Of  course,  for  poorer  efficiencies,  such  as  would  probably  occur  at 
the  lower  frequencies  considered,  the  data  rates  would  be  proportionately 
lower.  Although  we  consider  n/R-lo'^  to  be  a  reasonable  choice  for 
the  purposes  of  this  study,  other  values  could  equally  well  have  been 
used.  Fortunately,  the  main  conclusions  do  not  depend  strongly  on 
the  precise  value  of  n.  Moreover,  since  extensive  graphs  of  REP  are 


-54- 


glven  below,  total  power  requirements  can  be  obtained  for  any 
efficiency  and  data  rate:  simply  multiply  REP  by  R  and  divide  by  n. 

Figure  16  shows  REP  as  a  function  of  transmission  range  for  several 
frequencies  and  the  Keller  model.  For  comparison,  REP  for  the  Levin 
(dry)  model  Is  shown  for  a  frequency  of  1  Hz.  The  results  shown  In 
Fig.  16  apply  to  the  situation  where  the  receiving  antenna  Is  located 
near  the  center  of  the  waveguide — about  10-km  deep  for  the  Keller  model 
and  about  9-km  deep  for  the  Levin  (dry)  model.  Also  shown  In  Fig.  16 

(right-hand  axis)  Is  the  total  required  system  power  for  combinations 

*  -3 

of  efficiency,  n,  and  data  rate,  R,  such  that  n/R“10  .  Any  effects 

of  transmitter  burial  depths  are  included  Implicitly  In  the  efficiency, 

n. 

Several  important  conclusions  can  be  drawn  from  Fig.  16.  The 
higher  frequencies  considered  are  clearly  superior  to  the  lower  ones, 
as  evidenced  by  the  fact  that,  for  a  given  power,  far  greater  ranges 
are  achievable  at  5‘J  kHz  than  at  the  other  frequencies  shown.  This 
superiority  is  almost  entirely  due  to  the  quieter  noise  environment 
at  50  kHz  and  a  depth  of  10  km.  Reference  to  Fig.  13  shows  that,  for 
a  depth  of  10  km,  the  noise-transfer  function,  IY2I,  decreases  strongly 
and  monotonlcally  with  increasing  frequency.  Actually,  the  superiority 
of  the  higher  frequencies  Is  probably  even  more  clear-cut  than  shown  in 

_3 

Fig.  16 — it  should  be  easier  to  achieve  n/R“lO  at  50  kHz  than,  say,  at 
1  kHz,  since  antenna  efficiencies  typically  Increase  with  Increasing 
frequency. 

A 

We  do  not  Imply,  of  course,  that  n  will  be  Independent  of 
frequency  and  receiver  burial  depth. 


.  m. 
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For  purposes  of  discussion,  we  assume  that  10  megawatts  is  an 
acceptable  power  requirement  for  a  practical  system.  On  this  basis, 
transmission  ranges  on  the  order  of  1000  km  appear  achievable,  provided 
that  the  Keller  (or  Levin  (dry))  model  of  the  lithosphere  is  found  through 
field  measurements  to  be  a  reasonably  good  representation.  On  the  other 
hand,  the  system  power  requirement  based  on  the  Levin  (wet)  model  are 
so  high  for  transmission  ranges  of  interest  that  the  results  are  not  even 
worth  presenting. 

One  very  important  general  point  should  be  loted:  viz.,  the 
system  performance  depends  exponentially  on  signal-  and  noise-attenuation 
rates,  but  only  linearly  on  other  parameters.  Thus,  the  required  power 
at  a  given  frequency  depends  linearly  on  n  and  R,  and  exponentially  on 
the  transmission  range,  d.  Conversely,  achievable  transmission  range 
depends  only  logarithmically  on  efficiency,  data  rate,  and  power.  A 
misleading  impression  can  be  gained  if  syst^j.  performance  is  assessed 
on  the  basis  of  the  power  required  to  transmit  to  some  specified  range. 
Instead,  performance  should  be  assessed  on  the  basis  of  the  range  that 
can  be  attained  with  some  specified  power.  The  following  simple  numerical 
example  will  clarify  the  distinction;  As  indicated  above  fcr  the  Keller 
model,  a  range  of  1000  km  is  attainable  with  a  total  power  of  10  mega- 

_3 

watts,  provided  that  n/R*10  .  Assume,  for  purposes  of  discussion,  that 

future  analysis  shows  that  n/R=10  ^  is  the  best  that  can  be  achieved. 

If  the  1000-km  transmission  range  is  regarded  as  specified,  then  the  power 
requirement  increases  to  100  megawatts,  which  is  probably  impractical. 

On  the  other  hand,  if  the  10-megawatt  power  is  regarded  as  specified. 
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then  the  ten-fold  decrees.  In  efficiency* **  causes  the  e.erl.u.  tranenlsslon 
range  to  decrease  to  about  850  km-a  penalty  of  only  15  percent. 

Since  systee,  cost  ulll  depend  greatly  on  the  depth  to  which  antennas 
„ust  be  burled,  we  next  consider  the  penalty  In  transmission  ranjs 
(power  being  specified)  that  must  be  paid  If  receiving  antennas  are 
burled  shallower  than  the  nominal  center  of  the  waveguide.  Since 
practical  radio  coMnunlcatlon  systems  typically  have  many  more  receivers 
than  transmitters,  receiver  burial  depth  Is  probably  a  more  Important 

consideration  than  transmitter  burial  depth. 

Figure  17  shows,  tor  the  Keller  model  and  frequencies  of  10  and 
50  kHt.  the  dependence  of  required  power  on  transmission  range  for  several 
assumed  receiver  burial  depths.  For  a  frequency  of  50  kllr  and  a  total 
power  of  10  megawatts,  the  calculated  transmission  ranges  are  700-to-800 
km  for  receiver  depths  as  shallow  as  3  or  4  km.  These  ranges  are  only 
20-to-30  percent  less  than  those  corresponding  to  a  receiver  depth  of 
10  km.“  Since  drilling  costs  are  apt  to  be  a  large  traction  ot  the  total 
system  cost,  a  tradeoff  of  a  20-to-30  percent  reduction  In  transmission 
range  against  a  factor  ot  3  reduction  In  borehole  depth  would  seem 
advisable.  Such  a  tradeoff  cannot  be  carried  too  far,  however.  For  the 
Keller  model  and  a  frequency  ot  50  kllr.  use  of  a  receiver  de-.th  of  only 
2  km  would  cause  the  depth  degradation  function  (Figs.  U  and  15)  to 

*0r  equivalently,  a  ten-fold  increase  in  required  data  rate,  or 
ten-fold  decrease  in  allowable  system  power. 

**AEain  note  the  misleading  conclusions  that  can  be  drawn  if  trans- 
mission  lange.  rather  than  power.  Is  “"Fldered  to 

3  =rn:L  “ifaL;  L  increase  In  power 
of  nearly  3  orders  of  magnitude  (for  the  Keller  model). 
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becorae  so  large  that  the  corresponding  reduction  in  transmission  range 
would  probably  be  intolerable.  Note  also  that  for  lower  frequencies 
(e.g.,  10  kHz  or  less — see  Fig.  15),  the  fractional  degradation  in 
transmission  range  incurred  by  using  shallow  receiver  depths  is  less 
severe  than  at  50  kHz.  Nonetheless,  for  depths  greater  than  about 
three  kilometers,  the  higher  frequencies  are  still  preferable. 

Figure  18  is  completely  analogous  to  Fig.  17  except  that  it  applies 
to  the  Housley  BR  (assumed  transition  at  25  km)  model.  Should  this 
model  prove  to  be  a  reasonably  good  representation  of  crustal  conduc¬ 
tivities  over  large  geographic  regions,  then  transmission  ranges 
approaching  5000  km  could  be  achievable  for  receivers  placed  at  depths 
of  about  ten  kilometers.  Even  for  receiver  depths  as  shallow  as  3  or 
4  km,  ranges  in  excess  of  3000  km  are  a  possibility.  Again,  frequencies 
of  several  tens  of  kilohertz  are  preferable  to  lower  frequencies. 
Detailed  results  have  not  been  shown  for  the  Housley  EUS  model.  The 
achievable  transmission  ranges  based  on  this  model  would  be  of  the 
order  of  10,000  km. 

We  have  not  explicitly  considered  situations  for  which  the  litho¬ 
spheric  conductivity  profi.'e  exhibits  lateral  variations  over  the  propa¬ 
gation  path.  Of  course  such  variations  would  be  expected  to  occur  in 

situations  of  practical  Interest;  e.g.,  propagation  between  two  geo- 
* 

logical  provinces.  Mode  conversion  caused  by  lateral  variations  is 
not  expected  to  be  large  enough  to  alter  the  main  conclusion  of  this 
study  (Wait  and  SpieSj  2972).  Thus,  the  transmission  through  a 

A 

For  a  discussion  of  evidence  that  good  lateral  continuity  does 
exist  over  large  distances,  see  Heaaock  (1972),  pp.  1-9. 
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non-stratified  lithosphere  can  be  estimated  reasonably  well  by  simply 
summing  the  attenuation  suffered  in  each  region,  using  the  appropriate 
attenuation  rates  (Fig.  6,  p.  34)  and  path  lengths.  The  depth-degrada¬ 
tion  factor,  however,  depends  solely  on  the  properties  of  the  lithosphere 
in  the  vicinity  of  the  receiver. 

The  existence,  or  non-existence,  of  sharp  transitions  in  conduc¬ 
tivity  on  the  bottomside  of  the  waveguide  can  alter  the  maximum  trans¬ 
mission  range  by  about  a  factor  of  two.  An  example  of  this  effect  is 
given  in  Fig.  19,  which  shows  the  degradation  in  calculated  transmission 
range  that  occurs  as  the  depth  of  the  assumed  lower  reflecting  layer 
in  the  Housley  BR  model  is  Increased  beyond  25  km.  Of  course,  improve¬ 
ments  would  occur  for  shallower  layer  depths,  as  they  would  If  either 
the  Housley  EUS  or  SN  profiles  were  used  (Fig.  3,  p.  13).  Clearxy, 
field  experiments  (see  Hales,  1972;  Keller,  1972),  as  well  as  studies 
of  the  relation  of  laboratory  data  to  in  situ  rocks,  must  be  undertaken 
if  these  uncertainties  are  to  be  satisfactor i xy  resolved. 

In  regions  where  the  conductivity  profile  at  depths  shallower  than 
8  km  is  similar  to  that  shown  by  Keller  (Fig.  2,  p.  10)  rather  than 
as  shown  by  Levin  (Fig.  1,  p.  9),  antenna  burial  depths  of  only  3  or 
4  km  should  be  adequate  for  a  meaningful  propagation  experiment .  As 
Indicated  by  the  results  shown  in  Fig.  17,  other  conditions  would  also 
have  to  be  satisfied: 

1)  Combinations  of  transmitter  efficiency,  bandwidth  (R),  and 
total  power  should  be  selected  so  that  the  power  radiated 
into  the  lowest  TM  mode  is  at  least  a  few  watts.  (The 
corresponding  requirement  for  total  power  could  be  well  in 


excess  of  a  kilowatt.) 


Normalized  required  effective  power, — — ^  (Watts) 


Transmission  range  (km) 


Fig.  19— Transmission  Range  vs.  Effective  and  Total  Power  for 

Several  Depths  of  Assumed  Sharp  Transition  in  Conductivity 
(Housely  8R  Profile,  lO-km  Receiver  Depth) 


Total  required  power  for  '»7/R  =  10  (Watts) 
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2)  A  frequency  of  several  tens  of  kilohertz  should  be  used. 

As  discussed  by  Levin  (1971),  a  propagation  experiment  has  been  carried 
out  at  a  1  equency  of  2.2  kHz,  a  total  power  of  2.5  kilowatts,  and  with 
3-km  deep  boreholes.  The  negative  outcome  of  this  experiment  does  not 
provide  conclusive  evidence  that  favorable  conductivity  profiles  are  non¬ 
existent  in  nature.  The  2.2  kHz  frequency  was  much  lower  than  optimum, 
and  the  antenna  efficiency  might  well  h/.ve  been  too  poor  to  provide 
adequate  REP  from  the  2.5-kilowatt  total  power.  Further,  the  3-km  deep 
boreholes  may  simply  have  been  too  shallow  for  the  particular  region 
where  the  experiment  was  carried  out— deeper  boreholes  or  a  different 
test  site  being  needed.  Clearly,  very  careful  site  selection  and  antenna 
d-ssign  are  needed  for  a  conclusive  experiment. 
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This  report  presents  results  of  an  analysis  of  the  capabilities  that  might  be 
achieved  with  a  lithofspheric  communication  system.  The  analysis  uses  several  avail¬ 
able  models  of  contine'^^el  lithospheric  conductivity  as  inputs  to  calculations  of 
transmission  charactei  vtlcs  and  possible  communication  system  parameters.  Because 
of  prevailing  uncert*lncl«8»  significant  differences  exis.  among  the  various  models; 
and  a  wide  range  of  future  system  capabilities  is  found  to  lie  within  the  realm  of 

possibility. 

For  all  models  used,  the  c#»lculated  attenuation  rates  are  considerably  larger 
than  for  above-ground  transmission,  which  depends  on  propagation  in  the  earth- 
ionosphere  waveguide.  In  spite  of  this  fact,  reasonably  large  transmission  ranges 
might  be  possible  in  the  lithosphere  because  the  atmospheric  noise  fields  at  depth 
are  very  small,  having  suffered  heavy  attenuation  in  propagating  downward  from  the 
earth's  surface. 

The  major  factor  affecting  system  feasibility  is  whether  water-saturated  micro¬ 
cracks  exist  in  alJ  ‘^ajor  rock  types  at  depths  treater  than,  say,  5-to-8  km.  If  such 
microcracks  exist  assumed  in  one  model),  and  minimum  conductivities  in  the  crust 

exceed  10“5  -  10“^'  uihos/m,  then  the  calculations  show  clearly  that  a  practical  com¬ 
munication  system  is  not  feasible.  For  this  case,  the  attenuation  rates  are  so  high 
that  only  very  short  transmission  ranges  could  be  achieved,  even  if  large  amounts  of 
power  were  expended.  Field  experiments  will  be  needed  to  determine  with  certainty 
whether  saturated  microcracks  are,  in  fact,  present  at  the  depths  of  Interest. 
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However,  valid  reasons  exist  for  optimism  that  fluids  (or  fluid  films) 
should  not  be  an  important  factor  at  depths  greater  than  several 
kilometers,  and  that  conductivity  profiles  (at  depth)  based  on  labora¬ 
tory  data  for  dry  rocks  more  closely  approximate  actual  conditions. 
Further,  recent  laboratory  data  indicate  that  the  conductivity  of  dry 
rock  in  the  crust  could  be  much  lower  (10~8  to  10“^  mhos/m)  than 
previously  believed  (10"^  mhos/m).  It  is,  of  course,  difficult  to 
relate  data  obtained  from  laboratory  samples  to  conductivities  of 
Jock  in  situ.  Nonetheless,  calculations  using  conceptual  profiles 
based  on  the  assumption  of  dry  basement  rock  yield  results  that  are 
quite  encouraging.  For  example,  for  data  rates  of  a  few  tens  of  bits- 
per-second  and  power  expenditures  of  1-to-lO  megawatts,  transmission 
ranges  of  from  about  1000  to  many  tl ousands  of  kilometers  (depending 
on  the  profile  used)  appear  to  be  a  possibility. 

The  calculations  show  clearly  that  much  larger  transmission  ranges 
can  be  achieved  in  the  LF  band  (30-100  kHz)  than  at  the  lower  fre¬ 
quencies  considered.  This  ?ffect  is  due  mainly  to  the  heavy  attenua¬ 
tion  suffered  by  downward  propagating  atmospheric  noise  at  the  higher 
frequencies.  The  calculations  also  show  that  the  maximum  transmission 
ranges  are  somewhat  sensitive  to  the  details  of  the  model  profiles. 

For  example,  for  two  profiles  exhibiting  the  same  minimum  conductivity, 
the  transmission  range  (for  a  specified  power,  etc.)  can  vary  by  a 
factor  of  two  depending  on  the  conductivity  depth-gradients,  and  the 
precise  location  of  any  sharp  transition  layers  that  might  exist. 

The  depth  to  which  a  receiver  must  be  burled  to  achieve  satisfac¬ 
tory  performance  is  a  particularly  Important  parameter,  since  borehole 
drilling  costs  could  be  a  large  fraction  of  the  total  system  cost. 

The  transmiboion  ranges  given  above  apply  when  the  receiving  antenna 
is  located  near  the  "center"  of  the  waveguide — assumed  to  be  about 
ten  kilometers  deep  for  the  conductivity  profiles  used  in  this  report. 
For  profiles  representative  of  regions  stripped  cf  highly  conductive 
sedimentary  layers,  the  calculations  show  that  receiver  depths  as 
shallow  as  3  or  4  Km  could  be  used  at  only  a  15-to-20  percent  penalty 
in  transmission  range. 
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